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The thick brane model supported by a nonlinear spinor field is constructed. The different cases 

with the various values of the cosmological constant A ( = J are investigated. It is shown that 

regular analytical spinor thick brane solutions with asymptotically Minkowski (at A = 0) or anti-de 
Sitter spacetimes (at A < 0) do exist. 
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In the last years, the brane world models gained widespread acceptance for a description of the Universe (for a 
review, see [l|, Q). In such models it is assumed that our Universe is a four-dimensional hyperspace embedded in 
a higher-dimensional space. Such an approach allows to solve some important problems of a theory of elementary 
i— —i' particles. There are two types of brane models: thin and thick ones. In the first case, it is assumed that a thickness 
of the brane is infinitely small, and matter has delta-like distribution on the brane. However, from the physical point 
of view, it is obvious that this is just some kind of idealization, and more realistic models of the Universe demand 
introduction of the brane thickness. For this reason, a consideration of the thick brane models assumes ever greater 
, importance (see e.g. the recent review 0). 

Creation of models of the thick branes is quite interesting process by itself. The main difficulty consists in the 
1 necessity of obtaining a regular solution. Here we call the solution regular if the energy density per unit volume of 
the brane is finite. Usually, the thick brane solutions are being obtained as a result of the interaction of scalar fields 
with gravitational one. In this paper we present a new model of the thick brane supported by a nonlinear spinor field. 
It is well-known [J] that without gravitational field there exists a regular spherically symmetric solution for the 
, spinor field. It is reasonable to suppose that the inclusion of the gravitational field does not destroy these solutions. 
However, self-consistent asymptotically fiat solutions for a gravitating spinor field arc not known until now. Note 
that the solutions for a spinor field propagating on a curved background are known (for details, see Ref. Q). The 
cosmological solutions with a spinor field are known as well, see Refs. (6l|-[loj. In Refs. @-0 the cosmological 
solutions with the gravitating spinor fields are considered. In Ref. @ the classical and quantum evolution of a 
• • . universe in which the matter source is a massive Dirac spinor field and the universe is described by a Bianchi type I 
. £h ' metric is investigated. It is shown that there exists the transition from an Euclidean to a Lorentzian regimes. In Ref. 
(Toj a class of exact cosmological solutions with a neutral scalar field and Majorana fermion field was found. 
In this paper we consider a thick brane model supported by the nonlinear spinor field. The main purpose is to 
?3 ! show that in the model with the nonlinear gravitating spinor field the existence of localized regular multidimensional 
solutions trapping test matter fields is possible. 



II. 5D BRANE FROM NONLINEAR SPINOR FIELD 



We consider a five-dimensional gravitation with the nonlinear spinor field as a source of matter. Usually, in four- 
dimensional problems the following types of nonlinear terms are used: 
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Here we choose the simplest variant when the nonlinear term is taken in the form (ipip) ■ Then the Lagrangian of the 
spinor field will be 



C m = \(^V^\>- 0^) - rmptp + - (ipip) 2 ■ 
The corresponding five-dimensional Einstein and Dirac equations are 

R a A ~ \e a A R = *T a A + e /A, 
\iT a e a A D A -m + \(^)]^ = 0, 



(1) 

(2) 
(3) 



where a = 0,1,2,3,5 is the Lorentz index; A = 0,1,2,3,5 is the world index; e a is the 5-bein; r a are the 5D 
Dirac matrices in a fiat Minkowski space; Da4> = (pA — \u A ab Y a b) ip is the covariant derivative of the spinor ip; 
Tab = \ (raTb — r;T a ); V ^ = e al a ^ A 1 ^', m, X are some parameters; A is the cosmological constant. All definitions 
for the spinor differential geometry are taken from [ill ] . The energy- momentum tensor for the spinor field is 



rp A 



(4) 



\i> (T A e a B + T a g AB ) D B ^ + % -D B ^ (T A e a B + T a g AB ) j, + e A C m - 

l -^T A a Jli> + l -0T a A i>, 

where the last two terms vanish on-shell; T A — e A T a are the 5D Dirac matrices in a curved spacetime; g 
e a A e b B ri ah is the 5D contravariant metric tensor; rj ab = {+1, — 1, — 1, — 1, —1} is the contravariant metric tensor in the 

5D Minkowski spacetime; ip = ip^T® j s Dj rac conjugated spinor; Da^P — Tp (^d a + |w A ° & r a 6^ , where %p d a = dA^P- 

The 5D Dirac matrices in a flat Minkowski space are 



AB 



r u 

r 7 

p5 



hx2 
hx2 

-a 

a-, 



1 ,i = 1,2,3, 



-«I2X2 
il 2 x2 



(5) 
(6) 
(7) 



where 12x2 is 2 x 2 unity matrix, and are Pauli matrixes 

01 



1 

1 



-i 

1 



1 

-1 



We seek a wall-like solution for the system ©-Q. To do this let us choose the 5D bulk metric in the form 



ds 2 



<P 2 {r)ds 2 M -dr 2 



where ds 2 M is the 4D Minkowski metric. For the spinor field, we use the following ansatz 




(8) 



(9) 



Then, using Eqs. @ and ([3]), one can obtain 



->cXa 2 b 2 



A 

3' 



a' + 2 — a - ma + 2\a 2 b 

4> 

b' + 2—b + mb - 2Xab 2 



x 
3~ 

0, 
0. 



(a'b -ab') + \a 2 b 2 



A 

6' 



(10) 

(11) 
(12) 
(13) 
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Introducing new dimensionless variables f = mr, a — ay/2\/m, b = by/2\/m, x = x/(4A), A = A/m 2 , and using 
Eqs. (fl2|) and (fl3| for eliminating derivatives on the right hand side of Eq. (fTTj) . we have 



a + 2 z -a - a + a 2 h = 0, 



b' + 2—6 + b - &b 2 



2~~272 A 

r a b 3 - 



,/2 

6 2 3 V 



6 ' 



(14) 
(15) 
(16) 
(17) 



This system can be solved analytically. To do this let us find from Eqs. (fT6|) and (fT7|) the following first integral 

a6> 4 = C, (18) 

where C is an integration constant. Then Eq. (|15[) can be easily integrated giving three types of solutions: 
1) At A = 0. In this case 



1 
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t:C 1 + — kr 



a = ao exp < r 



b = b exp < -f + - 



1/4 



arctan 4 



f I + In I 1 + — x P 



16, 



— arctan 4W — f -In H x f 



(19) 
(20) 

(21) 



where ao, bo are integration constants. In the case under consideration do = &o = v C*. 
2) At A > 0. Here Eq. (JTSJ) gives 

1 1/4 



1 / / 2 4qA2+p2 

2 P± V P ^TT - 



(22) 



with 



,| = (,„, j : -i\l^(f + r ) 



4xC 



2xC 2 



P = 



A 



Q = 



A 



where fo is an integration constant. In this case, it is also possible to find analytical solutions for a and b but they 
are quite tedious, and that is why we do not show them here. 
3) At A < 0. In this case we have 



1 

4B 



[{B-p) 2 +Aq} 



1/4 



a = a exp if- -C, 



1 



— — arctan 

g|A| 



i 
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b = 6n exp < — f H — d — — arctan 

I 2 h\M 



2^9 
1 



2^9 



(B-p) 
(B-p) 



In 



-im 

2 



(ff-p) 2 + 4g 
4B 

{B-pf+Aq 
AB 



(23) 
(24) 

(25) 



with 



B = 7 exp 41 



4xC 



2kC 2 
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FIG. 1: The profiles of the energy density Tq for A = (the solid line), A = —0.5 (the dashed line), A = 0.5 (the dash-dotted 
line), x = 1. The profiles for the curves A = and A = 0.5 practically coincide. 



where ao,bo,^f are integration constants. 

The corresponding dimensionless energy density is 

in = — ^ln = —ixa o . 

One can see from Eqs. (fTS)) - (f2"4"|) that for A = and A < the energy integrals on the coordinate f are finite 

oo oo 

f°^^dr = -2xC 2 J ^<oo, (26) 

— oo 

i.e. the solutions are regular and one can expect that such solutions can be used for modeling the branes. It is just 
necessary to check up whether test fields will be trapped by these solutions (see a consideration of this question in the 
next section). On the other hand, in the case when A > this integral diverges, and that is why using such solution 
for modeling a brane is not straightforward (see discussion of this question in the Conclusion). 

To demonstrate the behavior of the solutions, let us solve the system p^|) - ([T7|) with the following boundary condi- 
tions 

1/2 



0(0) = 1, d>'(0) = 0, 5(0) = 6(0) = | 1 + y 1 - A ] . (27) 

The last condition follows from the constraint equation (fl~5|) . We set (f)(0) = 1 because we can redefine the Minkowski 
coordinates arbitrarily, <j)'(0) = because the bulk space should be symmetrical relative to the origin of coordinates 
r = 0. The spinor ([9]) had been chosen in the spinor representation, and in the standard representation the spinor 
components have to be either odd or even functions. That is why we choose 5(0) = b(0). The corresponding results 
for the energy density (|26[) for three different values of A are presented in Fig. [T] The behavior of the metric function 
<t> and the spinor functions a, b can be estimated from the analytical expressions obtained above. 



III. TRAPPING OF MATTER 



The five-dimensional localized wall-like solutions obtained above can be used for a description of a brane only if it 
will be possible to show that various test matter fields can be confined on such wall. As an example of such field, let 
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us consider here a test complex scalar filed x with the Lagrangian 

1 „ . „ „ 1 



L 



where tuq is the mass of the test field. Using this Lagrangian, we find the equation for the scalar field 



f 



c) 



dx_ 

dx B 



-mix- 



(28) 



Here \/— 5 g is the determinant of the 5D metric qab and x is a function of all coordinates x — x( xA )- Taking into 
account that the canonically conjugate momenta p^ = (E, p) are integrals of motion, we seek a solution in the form 

x(x A ) = X(f) expi-ip^). 
Substituting this ansatz in ([28]). one can find the following equation for X(f) 

X" + 4^X' + (p» Pfl -m 2 )X = 0, 

or, taking into account that p^p^ — <fi~ 2 (E 2 — ~^ 2 ), we have 

X" + d-X' + [(E 2 - f 2 ) cf>- 2 -m 2 ]X = 0, (29) 



where the prime denotes differentiation with respect to r. In the case when A = 0, the asymptotic behavior of Eq. ([T 
is 4> w f3x/\f\, where (3 = (§C^) 1/4 > 0. Then we have from ([29]) 



X" + -X' - mlX = 0. 

r 

This equation has an asymptotically decaying solution in the form 

p -m \f\ 
Jf(°) » D- . 



(30) 



where D is an integration constant and the index (0) refers to the solutions with A = 0. 

In the case when A < 0, the asymptotic follows from Eq. (|2"3")l : <\> w (-ji?) 1 ^ 4 . Then Eq. (|2^|) gives 

X" + 4\l ^-X' - m 2 X = 0, 

with an asymptotically decaying solution 

X^^Dexp 



(31) 



where D is an integration constant and the index (— ) refers to the solutions with A < 0. 

As a neces sary condition of trapping of the matter, one can require converging the field energy per unit 3- volume 
of the brane [12|], i.e. 



Etodx] 



J T»^r 9 dr= J 4 ^(E 2 + J? 2 )X 2 + m 2 X 

— oo — oo 

and also the norm of the field \ should be finite 

oo 

llxl| 2 = / V Z ^9X*Xdr= f°° ^X 



X' 



df < oo, 



(32) 



dr. 



Taking into account the asymptotic solutions (|30|) and (1311) for A = and A < 0, respectively, one can see that 
E to t and ||x|| converge asymptotically for both cases. Thus, it is obvious from the above analysis that the localized 
solutions which we found in section |TI] trap the test scalar field that indicates that such solutions may be interpreted 
as brane solutions. 
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IV. CONCLUSION 

In this paper, we have considered the ^-symmetric domain wall (thick brane) solutions supported by the nonlinear 
spinor field both with and without account of the five-dimensional cosmological A-term. It was shown that regular 
solutions, i.e. solutions with the finite energy density per unit 3- volume of the brane, do exist. In the case of positive 
A one can see that the total energy density over all space diverges. In the cases A = and A < the total energy 
density remains finite. It allows using such solutions for a description of the thick branes. To show this, we have 
considered trapping of the test scalar field on the wall. It was shown that asymptotically converging solutions exist 
for the test scalar field both for A = and A < cases. The asymptotic behavior of the brane solutions corresponds 
to the Minkowski (at A = 0) or anti-dc Sitter spacctimcs (at A < 0). 

As regards the case A > 0, the obtained periodic solutions for the energy density (see Fig. Q]) can be used for a 
description of the five-dimensional spacetime with the compactified fifth coordinate. In this case the size of the fifth 
dimension is defined by the period of the obtained solution for the spinor field. 
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